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Abstract 

An experiment demonstrating single-pixel single-arm complementary compressive microscopic ghost imaging based on a 
digital micromirror device (DMD) has been performed. To solve the difficulty of projecting speckles or modulated light patterns 
onto tiny biological objects, we instead focus the microscopic image onto the DMD. With this system, we have successfully 
obtained a magnified image of micron-sized objects illuminated by the microscope’s own incandescent lamp. The image quality 
of our scheme is more than an order of magnitude better than that obtained by conventional compressed sensing with the same 
total sampling rate, and moreover, the system is robust against intensity instabilities of the light source and may be used under 
very weak light conditions. Since only one reflection direction of the DMD is used, the other reflection arm is left open for 
future infrared light sampling. This represents a big step forward toward the practical application of compressive microscopic 
ghost imaging in the biological and material science fields. 
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I. INTRODUCTION 


Ghost imaging (GI) emerged two decades ago as a technique based on the correlated nature of both classical and 
quantum fields. Initially, ghost images were formed from two correlated light fields and two photodetectors: a high 
spatial-resolution detector used to record one light field which had not interacted with an object to be imaged, and 
a bucket detector without spatial resolution which was used to collect the light coming from the object. Neither 
detector can produce an image on its own, but a “ghost” image can be constructed by correlating the signals of these 
two detectors. The first demonstration of ghost imaging (GI) El! utilized biphoton pairs which were produced by 
spontaneous parametric down conversion in a nonlinear crystal, hence the phenomenon was interpreted as the result 
of quantum entanglement of the photon pairs [31. Subsequently, further theoretical and experimental work showed 
that GI is also achievable with pseudo-thermal [4-[T0j as well as true thermal light [Tl|, 13 an d can be explained with 
a classical statistical model EMI. 

Later, it was realized that the purpose of the reference beam was merely to measure the field distribution at the 
object, so it could just as well be replaced by a single predetermined spatially modulated object beam plus the bucket 
(single-pixel) detector. This modified technique was called computational GI by Shapiro [15[ and demonstrated by 
Bromberg et al. (16). But a few months earlier, the same optical protocol had been proposed and demonstrated 
independently by Baraniuk et al., who combined it with the compressive sampling (compressed sensing) algorithm 
E3-E3 , and called it “Single-pixel imaging via compressive sampling” [13, Hl| . The only difference between single¬ 
pixel compressed sensing (CS) imaging and computational GI was that the spatial modulator was placed behind the 
object instead of in front. 

It is known that traditional GI reconstruction algorithms are linear while CS algorithms are nonlinear. In CS, 
optimization is applied to recover a signal from incomplete or noisy observations of the original signal through random 
projections, provided that the signal is sparse or compressive in some basis. The application of CS in the field of 
imaging allows one to retrieve high resolution images from many fewer measurements than those established by the 
Nyquist-Shannon criterion. Recently, the GI community has employed CS to obtain compressed ghost images [22j| . 
greatly reducing both the acquisition time and the number of frames required. These improvements have motivated an 
ongoing effort to implement technologies based on GI such as multispectral imaging (2j] , optical encryption E3, [13 , 
and most recently, adaptive compressive GI (26) . Although GI and single-pixel compressive imaging were historically 
independent and developed in parallel with each other, both techniques promise a resource-efficient alternative to 
array detectors, permitting us to reduce operational problems involved in systems based on raster scanning. However, 
practical random patterns sequentially fed into a DMD only have two values, 0 or 1, corresponding to the ±12° angles 
at which the micromirrors are deflected, which is not the best range of values for compressive GI. It is interesting to 
note that when both reflected beams from the DMD are used to reconstruct a ghost image, the corresponding frames 
are complementary. Accordingly, instead of wasting the 0 (—12° reflected) beam we have developed a novel technique 
which we call complementary compressive imaging, that makes full use of both reflections to dramatically improve 
the image quality [27l |. 

On another note, some time ago Luo et al. mm reported a technique that they called correspondence imaging, 
in which a positive or negative image is retrieved by conditional averaging of the reference signals; that is, only those 
reference data that correspond to positive or negative intensity fluctuations of the bucket signal are selected for simple 
averaging, without the need to multiply by the bucket detector intensity itself. It seemed to defy intuition in that no 
direct second-order correlation calculation is performed, while compared with conventional GI the number of exposures 
used to reconstruct the images and consequently the computation time are greatly_jeduced. Further developments 


of this positive-negative image concept have been presented in many other papers 


]. By alternately measuring 


the bucket values of a random binary pattern and its inverse, Sun et al. |33|, (34( reconstructed a ghost image by 
correlating non-inverted patterns with the differential signals of the complementary illuminated pairs, and normalized 
the bucket signal with respect to the positive or negative intensity fluctuations averaged to 0. Their method was 
essentially the same as that of correspondence imaging, but their GI bucket values were used as a series of weighting 
factors instead of just 0 and 1. Their projected inverse patterns were only used in generating the differential bucket 
values and not in the correlation computation, so their scheme is not suitable for compressive GI, which is a nonlinear 
approach and needs a one-to-one relationship between the modulation patterns and measured values. Moreover, it is 
very difficult to project black and white speckles onto tiny biological tissues by a projector. To our knowledge, the only 
previous report on the use of GI in actual microscopy was based on entangled photon pairs, which is a complicated 
and expensive light source |35j |. Ghost imaging with a classical thermal source was studied in a two-arm microscope 
imaging system but it was only a theoretical simulation with a simple double-slit object [36j. Recently, CS was also 
applied to microscopy (371, 38); in (38j, Studer et al. tested their system on a sample of fluorescent beads which were 
sparsely distributed. In addition, they used binary patterns of a shifted and rescaled form (1 + h)/2 where h is a 
Hadamard sequence, thus each entry of patterns is either 0 or 1, still not the best range of values. In (39), Radwell et 
al. proposed an “adaptive” mode that provides the benefits of both high frame rate and high resolution, but still used 
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Hadamard patterns for image reconstruction of a simple silicon CMOS chip. Here, as an alternative to complementary 
measurements in both reflection arms as in j2?| , a prototype single-pixel microscope with complementary modulation 
is presented. It can not only image complex gray samples but also produce very satisfactory image quality, more than 
an order of magnitude better than traditional compressive imaging. Moreover, since only one reflection direction of 
the DMD is used, the other reflection arm is left open for future infrared light simultaneous sampling (which has just 
been demonstrated by Radwell et al. f33|) . 

On the basis of our former work [27H29J, we have extended our complementary compressive imaging technique to 
the imaging of real biological samples using an ordinary microscope’s own illumination lamp. In our proof-of-principle 
experiment, the object’s image is projected onto a DMD, which is encoded with randomly modulated array patterns 
alternated with their inverse (complementary) patterns. Moreover, we use only one reflection arm with a single-pixel 
ultra-weak light detector rather than the two detectors in [27j], leaving room for future simultaneous infrared light 
detection in the other arm. Through our method, complex gray-scale images of the mouthparts of a female mosquito 
have been reconstructed with far fewer measurements and a performance more than an order of magnitude better 
than that of traditional CS, opening exciting prospects for the important field of microscopic imaging. 


II. EXPERIMENTAL SETUP AND IMPLEMENTATION 

The experimental apparatus is illustrated in Fig. [l] The thermal light from a halogen lamp illuminates the sample 
table from overhead with an ordinary incandescent lamp (a noncoherent thermal source). Light from the object 
passes vertically down through the objective, and after reaching a flippable beam splitter may be viewed directly 
through the eyepiece or transmitted by various mirrors and lenses to the DMD. A 10x magnification is achievable 
with just the objective, which has a numerical aperture of 0.25. The DMD consists of 768 x 1, 024 micro-mirrors each 
of size 13.68 x 13.68 /irn 2 . Each mirror rotates about a hinge and can be switched between two positions oriented at 
+12° or —12° with respect to the plane of the DMD; the +12° position corresponds to bright pixels 1, and inversely, 
— 12° to dark pixels 0. The modulation frequency of the DMD can reach 32.5 kHz. A photomultiplier tube (PMT) 
(Hamamatsu H7468-20) is used as the bucket (single-pixel) detector for collecting the total light intensity reflected 
from the DMD, with some neutral density filters in front for suitable attenuation. Therefore, our proof-of-principle 
single-pixel microscopic imaging experiment is performed under ultra-weak light conditions where ordinary microscopy 
photography would be impossible. Generally, in biological microscopy, the object is a fluorescing sample like cells, 
tissues and beads, the light from which is extremely weak. Our sample consisted of the mouthparts of a female 
mosquito, of about 0.8 x 2.4 mm 2 in size. In order to show what the object is and to make sure the image is clear, 
we first took a photograph of the mosquito sample (as illustrated in the upper left corner of Fig. []} as a reference 
image directly in front of the eyepiece of the microscope with an ordinary charge coupled device (CCD) monitor and 
a common flippable beam splitter, under bright illumination. 



FIG. 1. Schematic of the experimental setup. 

In our configuration the size of the image to be reconstructed is 768 x 1,024 pixels, equal to that of the DMD. 
However, the time consumed by CS increases exponentially with the dimensions, so to recover such a large image 
we performed a row-scanning strategy: the DMD was activated row by row and the same number of complementary 
exposures were taken for each row. Since the light source can be considered to be stable within two reasonably short 
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modulation time intervals, we modulated the complementary frame pairs of each row alternately. This was repeated 
until all rows had been measured. Then we calculated the differential values of adjacent elements (in a one-to-one 
complementary relationship) of the bucket signal with odd and even index numbers, corresponding to complementary 
differential frames. After all rows had been reconstructed, the whole image was retrieved via the image mosaic method. 


III. COMPLEMENTARY COMPRESSIVE GHOST IMAGING FRAMEWORK AND RESULTS 

In conventional GI, an image of pixel size p x q can be acquired by correlating the reference frame intensities Ir 
with the bucket detector signal Sb- 

GI = (S b Ir(x r )), (1) 

where (• • •) is the average operator, and xr denotes the spatial position in the reference detector frames. Each bucket 
detector value can be defined as 

Sb= Ir(xi, Xj)T(xi, Xj)dxidxj, (2) 

J o 

where Xi, x 3 denote the row and column indices, respectively, T represents the intensity transmission function of the 
object, and O is the area of the beam. 

In compressive GI, we can flatten the two-dimensional image into a one-dimensional column vector, and denote it 
by x of length N, where N = p x q. Each pattern of the DMD can also be reshaped into a row vector, and then M 
such vectors re-arranged into a measurement matrix A of size M x N. Now the sequence of bucket values is also given 
as a column vector y of size M x 1, thus the problem can be described by the set of linear equations: 

y = Ax + e, (3) 

where e of size M x 1 denotes the noise. When M > TV, these equations can be easily solved by various iterative 
methods, such as lower upper decomposition, Gauss-Jordan elimination, Cholesky decomposition, or singular value 
decomposition. However, as the image size increases, the computational overhead grows rapidly. It is known that 
a natural image generally has a sparse representation in a certain basis T (e.g. Haar wavelet basis, discrete cosine 
transform basis, or Fourier transform basis). This prior knowledge is very helpful for addressing an ill-conditioned 
problem where the number of measurements is fewer than the total number of pixels N, and for reducing the number 
of measurements required for the image acquisition. Here we expand x in an orthogonal basis ’k = [ipi,tp 2 ,..., ipN~] as 
follows: 

N 

x = Ta/, or x = ^^x'ii/ji, (4) 

i —1 

where x 1 of size N x 1 is the coefficient sequence of x. We say that x' is fc-sparse if |{i : / 0,1 < * < N}\ < k. 

Then the problem becomes 

y = A^x' + e. (5) 

For faithful image reconstruction, rows of A should be incoherent with the sparse basis. It is found that a completely 
random matrix A works well, from which accurate recovery is possible by convex optimization fl9j |: 

min i ||y - A$x'\\l + r||a; , || 1 , (6) 

where r is a constant scalar weighting the relative strength of the two terms, and ||- ■ || ( , stands for l p norm defined 
as (|M| p ) P = Eti N P - The first term is small when the optimal x' is consistent with Eq. © within a small 

error. The second term is small when x' is sparse. Recent research (40j has proved that the use of total variation 
(TV) regularization instead of the l± term in CS problems gives a sharper recovered image by preserving the edges or 
boundaries more accurately, and the gradient of an image is generally sparse as well. Therefore, we use the TVAL3 
solver [ic| here. 

According to the method proposed by Sun et al. 33. H3| , ghost images can be retrieved by correlating differential 
signals A Sb and non-inverted patterns Ir(xr). That is, GIc = (ASbIr(xr)). The differential bucket signal in the 
second-order correlation form can be viewed as weighting factors, which indirectly gives the correlation coefficient 
between Ir(xr) and the object (i.e., A Sb > 0 and A Sb < 0 indicate a positive and a negative linear relationship, 
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respectively, and A Sb = 0 means no correlation), thus their method has a performance somewhat better than 
traditional second-order correlation but actually much poorer than that of compressive GI with the same number of 
measurements. Although the problem to be solved whether using intensity correlation or CS nonlinear algorithms is 
the same, the former correlation fails to work in CS because it destroys a precise correspondence relationship. 

We now consider the complementary sampling scheme proposed in our previous work [271 ]. Since each complementary 
frame pair appears alternately, the odd and the even subscripts constitute two supporting mutually related sets: 
odd = {1, 3,5,...}, and even = {2,4,6,...}. We denote the measured bucket signal and the measurement matrix used 
as y' and A !, respectively, then the complementary differential bucket signal is defined as 


V 2 (y°dd Veven)' (^) 

Here, y decreases the variance of the independent and identically distributed noise by a half [27). Accordingly, the 
complementary differential frames can be written as 


A — ~ ( A'odd - A' even ) ~ ^odd ~ 9 > 

Z ZMxN 


( 8 ) 


where the inverse frames A’ even = 1 mxn — A ! odd , % MxN stand for an array of all and we have let A! odd and A' even 
denote the submatrix of A 1 that contains all the columns corresponding to the support sets odd and even. Notice that 
A a and A e on the DMD are binary matrices consisting of the two values 0 or 1, and thus A is a ±0.5 binary matrix, 
which generates the “positive-negative” light modulation. Then Eq. m can be rewritten as 

liy'odd-y'even) = ( A 'odd~i )^' + \ ( e l ~ e 2 )■ (9) 

2 2MxN 2 


By utilising this technique, we have reconstructed the image of the mouthparts of a female mosquito, as pictured 
in Fig. [Ha), in which even the flagellum, clavola, upper and lower jaws, and labium of the mouthparts can be seen 
clearly. For comparison, we also show the image recovered via traditional CS (see Fig. [2])b) ) with the same total 
sampling rate of 50 % (that is, complementary CS samples 256 pattern pairs (i.e., 512 patterns in total) for each row 
of length 1024, compared to 512 patterns for conventional CS). From the experimental results, we can see that the 
image quality of our new method is much better than that of conventional CS. 



FIG. 2. Reconstructed images of the mosquito mouthparts of pixel size 768 x 1,024. Image retrieved by (a) complementary 
CS, and (b) conventional CS. 
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IV. PERFORMANCE AND DISCUSSION 


In our earlier work ]27j], we already proved by simulation that the total number of measurements required in the 
two-detector complementary sampling scheme is much smaller than that in traditional CS methods for the same mean 
square error (MSE) value, although the former needs to measure the signals in both reflection arms of the DMD. 
Moreover, it was also demonstrated that the computation time is shorter than that of traditional CS. As the system 
proposed in this paper is an extension of our previous scheme to the microscopic regime but with only one detector 
arm, it embodies all the same advantages. Since the modulation frequency used here is 450 Hz and the total number 
of measurements is the same, the total acquisition time for both complementary CS and conventional CS is altogether 
14.56 minutes. However, for the same image quality, conventional CS will require a significant increase in the number of 
exposures and consequently measurement time. Furthermore, an added bonus is the robustness of complementary CS 
against any instability of the light output from the lamp of an actual microscopic imaging system, which would affect 
the bucket detector signal y'. To simulate the deleterious effect of unstable illumination, we programmed the DMD 
such that the central row (Row 384) displayed 616 random patterns repetitively, as a result of which the measured 
bucket detector signal showed corresponding periodic variations (see Fig. [Ha)). This instability would result in poor 
image quality with conventional CS. However, by the pairwise subtraction algorithm within each complementary pair, 
the bucket detector values are converted into a “positive-negative” intensity distribution vector with a mean of ~ 0, 
as shown in Fig. [3|b), thus compensating for any undesirable fluctuations of the average illumination intensity. 



FIG. 3. Simulation of the effect of intensity instabilities on the bucket detector fluctuations: (a) y' and (b) y. 

To obtain a more quantitative measure of the reconstructed image quality, the peak signal-to-noise ratio (PSNR) 
is defined here as a figure of merit: 

PSNR = 10 log(255 2 /MSE), (10) 

where 

MSE = - ^ , [U 0 (r, c) - U(r, c)] 2 . (11) 

pq zJ r,c=l 

This MSE describes the distance squared between the reconstructed image U and the original image U 0 for all p x q 
pixels. Naturally, the larger the PSNR value is, the better the quality of the image recovered. 

In Fig.[]]the PSNR values of complementary CS and conventional CS are plotted as a function of the total sampling 
rate, which is defined as the ratio between the total number of measurements and the length of the image signal. 
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FIG. 4. PSNR values of complementary CS and conventional CS versus the total sampling rate. 

the 
the 
the 

same total sampling rate, i.e., the number of measurements is significantly reduced for the same image quality. The 
maximum difference in the PSNR is more than 15 dB, i.e., corresponding to a maximum MSE difference of about 32. 

Since the values of A a and A e all consist of either 0 or 1, and obey the Bernoulli distribution, the mean of both 
matrices is always close to 0.5. The impact of the background signal yb = \ MxN ^x' (which can also be viewed as 
the background noise) on image quality should not be neglected. Indeed, our Eq. © subtracts this background term 
and gives a better performance than Eq. ©■ Moreover, the mean of A is close to 0, which in fact helps to pick out 
useful information from the bucket signal. 


As shown in the figure, (yellow circles: complementary CS, pink diamonds: conventional CS), for both methods 
image quality increases with the sampling rate, tending to saturation at a high level. It is clearly visible that 
image quality of complementary CS is more than an order of magnitude better than that of conventional CS with 



FIG. 5. Schematic diagram of the DMD. The micromirrors in the +12° (—12°) orientation are shown in white (black). The 
lilac and peach colored areas represent the modulated and the nonactive regions of the DMD, respectively. 


Another source of noise which we should consider is the DMD. As illustrated in Fig. [5J the random patterns are 
only encoded onto the activated area a of the DMD, while all the other micromirrors in the remaining region b are 
in the —12° position. The gaps between adjacent micromirrors and gaps incurred during flipping of the micromirrors 
which expose the metal surface beneath will incur some random noise, including mirror-like and diffuse reflections. 
The former will not be reflected into the direction of detection and so can be neglected. However, the noise from dark 
pixels in area b and from gaps generated by micromirror flipping within the image area is worthy of attention. To 
quantify this noise, we need to add the following extra eight measurements: 
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where a, b denote the signals from modulated and the nonactive regions of the DMD, and the subscripts “+” and 
” denote the states when all micromirrors are aligned in the +12° or —12° positions, respectively. The superscript 
“ object ” means that the object lies within the detector beam path, otherwise it is outside. Here /i,/ 2 , • ■ •,/s are 
the corresponding measured values of the eight equations. By solving these equations, we can easily calculate b ob ^ ect , 
c°bject, an d a°y ect . Considering the dark noise Cdark of the PMT, we can rewrite Eq. m as 

object 

y = AVx 1 + (c dark + b°y ect + 0.5—- a°y ect ), (13) 

a+ 

^object 

where t = — h —— denotes the total intensity transmission of the object. Sinse the DMD pixels are randomly modulated, 
only half of the light flux in area a is detected in each exposure, so we only need discuss the noise in half of the (object’s) 

■ q object 

transparent part. Except for c dark which is a variable, the other two terms (b°_ jec and 0.5—— a°_ Jec ) in the noise 
e are all constants, and thus can be subtracted. However, even after these two items are subtracted, the background 
noise i A/xAr 'kx / mentioned above still plays a primary role in image degradation and the image quality is not expected 
to be improved dramatically. Since the dark noise is independently and identically distributed in both y' odd and y ' even , 
complementary CS can not only decrease the variance of the dark noise by half but also remove the impact of the other 

two kinds of noise (6_ ject and 0.5via Eq. ©. Although the complementary differential bucket signals 
seem to be a subtracted and rescaled version of those in conventional CS imaging, our method actually averages the 
variance of the noise, thus significantly improving the SNR in the measurement process. 


V. CONCLUSIONS AND PROSPECTS 

In summary, we have proposed and experimentally demonstrated a single-pixel single-arm DMD-based compressive 
ghost microscopy system, based on direct complementary modulation of the DMD. In order to incorporate the micro¬ 
scope into this system, the image is directly focused onto the DMD plane, rather than using structured illumination by 
a digital projector or two-arm sampling of traditional GI. The mouthparts of a female mosquito have been successfully 
imaged and identified. An image quality more than an order of magnitude better than that obtained by conventional 
CS with the same total sampling rate has been realized, and moreover, the system is very robust against intensity 
instabilities of the illumination source, while imaging is still achievable under very weak light conditions. Compared 
with the ±12 ° double-arm detection scheme in our previous work |2?| . it allows room for future simultaneous near- 
infrared light detection in the other reflection arm. Our scheme should open exciting prospects for microscopy in the 
biosciences and material sciences, notably for biological imaging of sparse fluorescent objects such as cells and tissues. 

Note: 

Here it may be interesting to recall the submission history of our paper. We first submitted our manuscript to 
Appl. Phys. Lett. (Doc. ID: L14-04737) on 25 May 2014. The editors told us that it did not meet the timeliness 
requirement for rapid publication in their journal, so we submitted the manuscript to Opt. Express (Doc. ID: 213211) 
on 31 May 2014. We received the first review decision on 30 June 2014, which was rejection on the grounds that the 
earliest idea of using a DMD instead of an SLM in computational ghost imaging was reported in the paper by Sun 
et al. [33;], and the same method of using opposite light pattern pairs in another conference paper by Sun et al. [34|, 
where they acquired the differential signals of each illumination pair for correlation calculation, rather than absolute 
signals. Actually, the use of a DMD instead of an SLM in computational GI is nothing new and was reported in many 
papers long before Sun’s, e.g. (23l - [25j . In addition, as mentioned in Introduction, Sun’s method was essentially the 
same as that of “correspondence imaging”, first proposed by Luo et al. in [28], [29j , but their GI bucket values were 
used as a series of weighting factors instead of just 0 and 1. 

We then modified the content of our manuscript and submitted it to Optics Letters (Doc. ID: 216782) on 10 
July 2014. We received the first manuscript decision on 25 July 2014. This time, a referee recommended that we 
should explain more clearly how and why the complementary pattern strategy is better than prior work on the linear 
reconstruction of images like in traditional GI, or the nonlinear approach of single-pixel imaging. 

On this advice, we added substantially more content to our manuscript, which then greatly exceeded the length 
restrictions of Optics Letters, so we once again submitted it to Optics Express (Doc. ID: 224607) on 7 Oct 2014. 
We received the first decision on 28 Oct 2014, but strongly disagreed with the reviewer’s comments and appealed 
for another review, which was granted on 5 Nov 2014. We received the review decision on 7 Jan 2015. The new 
reviewer suggest that the manuscript can be published only after we clear up the picture and the presentation of 
“ghost imaging”, and present it as the compressive imaging aspect only. However, we think that the only difference 
between single-pixel compressive imaging [20] and computational GI jl5[ was that the spatial modulator was placed 






behind the object instead of in front. Furthermore, the GI community has already employed CS to obtain compressed 
ghost images (see in i22,;23y25:j,:26j, Opt. Lett. 37 , 1067-1069, 2012, Phys. Lett. A 376 , 1519-1522, 2012, Phys. Lett. 
A, 377 , 1844-1847, 2013). Most intriguingly, the CS community also developed single-pixel compressive architecture 
based on an active illumination concept, as demonstrated in [37| and Appl. Opt. 405 , 405-414, 2011, where the 
spatial modulator was placed in front of the object, just as the same with the optical layout of computational GI. 
Removing the compressive aspect, i.e. taking the measurement matrix to be full rank, this imaging scheme turns into 
a mechanism named multiplexed illumination (in: Proc. IEEE International Conference on Computer Vision, 2003, 
808-815). For example, in Hadamard-nmltiplexed illumination, about half of the light sources operate simultaneously, 
creating brighter, clearer captured images. In fact, the mathematical models of single-pixel compressive imaging, 
computational GI or multiplexed illumination are the same [22j], only with the change of reconstruction algorithms, 
which will in turn result in the difference of the number of measurements. Although GI, single-pixel compressive 
imaging and multiplexed illumination were historically independent and developed in parallel with each other, these 
techniques promise a resource-efficient alternative to array detectors, permitting us to reduce operational problems 
involved in systems based on raster scanning. 

Besides, Padgett’s group should be congratulated on quickly realizing the dual infrared and visible microscope 
system, as described in |39j . 
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